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Defense  Research  Laboratory'  is  the  designation  used  in 
this  report.  On  1 October  1968  the  title  was  changed  to 
Applied  Research  Laboratories. 
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(U- !*0W)  The  efficiency  in  computing  the  Cepstrum  of  narrow  band 

signals  has  been  increased  by  the  use  of  direct  quadrature  sampling 
of  narrow  band  signals.  As  a practical  application  of  Cepstrum  to 
sonar,  this  is  an  important  step.  It  means  that  a Cepstrum  can  be 
computed  rapidly  on  a general  purpose  computer  without  the  use  of 
fast  special  purpose  hardware.  As  soon  as  the  programming  is  per- 
fected, the  new  method  will  be  used  to  analyze  more  real  data.  One 
of  the  purposes  will  be  to  test  the  ping-to-ping  stability  of 
selected  submarine  data. 


r 


( U- P0H9 ) More  work  was  done  on  a vocoder  analysis  of  active  sonar 

returns.  This  work  lias  been  interesting,  but  due  to  a lack  of  time 
it  will  have  to  be  discontinued. 


(U  A very  efficient  program  to  compute  quadrature  cross correlations 

based  on  the  fast  Fourier  transform  ( FFI)  was  developed.  A figure 
is  presented  showing  its  computational  speed  versus  the  number  of 
lags  computed  (see  Fig.  7). 
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II.  CENTRUM  ANALYSIS 
(J.  A.  Shooter) 

(C)  In  Quarterly  Progress  Report  No.  8 on  Contract  NObsr-95181  (u), 

it  was  shown  that  the  power  spectrum  of  a narrow  band  signal  can  be 
computed  directly  from  the  signal's  quadrature  components.  The 
computation  time  thus  saved  was  approximately  a factor  of  l6.  This 
technique  has  now  been  successfully  applied  to  computing  the  Cepstrum 
of  a narrow  band  signal.  An  illustration  is  shown  in  Fig.  1,  the 
data  in  this  case  being  a return  from  a stern  aspect  submarine.  The 
transmitted  signal  was  a linear  FM  sweep  with  a bandwidth  W of  625  Hz, 
a time  bandwidth  product  TW  of  64,  and  a carrier  frequency,  f , of 
5000  Hz.  Each  quadrature  component  was  sampled  at  a rate  of  2W.  The 
power  spectrum  of  the  positive  frequencies  is  shown  in  Fig.  la.  The 
frequency  axis  begins  at  f - f , and  f increases  to  f = f^  + W at 
the  center  of  the  plot.  Then  f changes  to  f - f^  - W and  increases 
to  f^  at  the  end  of  the  frequency  axis.  This  rather  strange  frequency 
axis  is  in  keeping  with  the  Cooley- Tukey  version  of  the  fast  Fourier 
transform,  and  it  is  necessary  to  leave  it  in  that  order  so  that  the 
next  Fourier  transformation  will  put  the  Cepstrum  peaks  in  the  correct 
position.  This  is  not  a symmetrical  function. 

( U- P0UO ) The  next  step  in  computing  a Cepstrum  is  shown  in  Fig.  lb,  which 
is  a filtered  logarithm  of  Fig.  la.  In  this  case,  it  was  assumed 
that  the  log  power  was  a time  series  sampled  at  20  kHz.  The  filter 
was  a recursive  band-pass  1-pole  Butterworth  filter  set  at  550  to 
6000  Hz.  The  filter  frequency  limits  are  related  to  the  time  delays 
between  superimposed  waveforms  or  target  size.  The  upper  and  lower 
limits  cover  target  sizes  from  stern  to  beam  aspects.  In  addition, 
to  remove  ringing  due  to  the  filter,  the  first  and  last  10^  of  the 
filtered  log  power  was  smoothed  to  zero  by  half  cosine  bells. 
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(c)  QUADRATURE  CEPSTRUM 


FIGURE  1 

COMPUTATION  OF  A CEPSTRUM  IN  QUADRATURE 
FOR  STERN  ASPECT  SUBMARINE  DATA  (U' 
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(C)  The  last  step  shown  in  Fig.  lc  is  the  square  of  the  modulus  of 

the  Fourier  transform  of  Fig.  lb.  In  addition,  a Hanning  window  was 
applied  to  the  real  and  imaginary  components  before  squaring.  This 
final  result  is  the  Cepstrum.  For  512  quadrature  components,  each 
sampled  at  twice  the  bandwidth,  the  computation  time  for  the  Cepstrum 
was  about  3 sec.  This  was  using  FORTRAN  language  on  a Control  Data 
Corporation  3200  computer.  In  machine  language  on  an  ordinary 
general  purpose  computer  there  should  be  no  trouble  in  computing  a 
Cepstrum  in  less  than  1 sec . The  quadrature  Cepstrum  is  actually 
the  envelope  of  the  Cepstrum,  and  the  results  are  interpreted  in 
the  following  manner.  The  main  peak  is  a correlation  of  a strong 
highlight  midway  between  the  stern  and  the  sail  with  the  highlight 
from  the  sail.  The  correlations  between  the  highlights  of  stern  and 
of  the  bow  are  very  weak,  and  they  are  almost  lost  in  the  noise . 

( U- SQOfti)  These  results  are  quite  encouraging.  The  next  few  experiments 

will  determine  the  proper  data  lengths  to  prevent  "wrap  around"  and 
to  optimize  the  signal- to- noise  ratios.  Once  the  data  lengths  are 
settled,  then  a ping-to-ping  stability  test  will  be  carried  out  to 
determine  the  suitability  of  Cepstrum  as  an  aid  in  classification 
of  submarine  echoes . 


(The  reverse  of  this  page  is  blank. 
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III.  VlWPKR  ANALYSIS  OF  ACTIVE  SONAR  DATA 
(J.  A.  Shoot  or  mid  R.  E.  Sent erf itt) 

[n  ^PK  No.  8 on  Contract  NObsr-98181  (u),  an  effort  was  made  to 

display  fin  TO  return  on  a time  vs  frequency  plot.  The  results  fit  that 

time  were  pur.? line*.  It  Is  now  known  that  the  difficulty  was  in  the 
fact  that  the  overlapping  TO  waveforms  can  Interfere  to  produce  peaks 
In  the  plots.  This  is  demonstrated  In  Fig.  ?,  which  shows  a vocoder 
plot  of  four  closely  spaced  overlapping  waveforms.  In  these,  art i- 
flclal  data  are  used  and  each  of  the  waveforms  is  identical.  It  is 

seen  that  the  output  of  each  filter  Is  unique,  and  therefore  the  peaks 

should  not  he  Interpreted  as  highlight  structure. 

(H-IW'  This  is  not  to  say  that  a time- frequency  plot  could  not  be 

useful.  Figures  *<,  4,  p , and  o show  vocoder  plots  of  a bow  aspect 
submarine.  These  represent  four  consecutive  pings,  and  the  results 
ft. tv  consistent.  At  this  time,  the  peaks  have  no  meaning,  other  than 
to  any  that  the  target  is  complex.  The  high  frequency  ripple  super- 
imposed on  the  large  peaks  are  caused  by  frequency  leakage  through 
a rather  poor  filter,  and  they  have  no  meanings. 

{ This  could  be  an  Interesting  problem  as  related  to  echo 
classification,  but  time  does  not  allow  its  perusal. 


UNCLASSIFIED 

— "f>»  "9Min>l  Only 


BAND-PASS  CENTER  FREQUENCY  - kHz 


UNCLASSIFIED 


START  OF  FM  SLIDES 


5 359 


A 

rrT~r 


FIGURE  2 

VOCODER  ANALYSIS  OF  ARTIFICIAL  DATA 
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FIGURE  3 

VOCODER  ANALYSIS  OF 
BOW  ASPECT  SUBMARINE  DATA  (U) 
TW  « 64  W = 625  Hi  fQ  * 5000  Hi 
SEQUENCE  13  TAPE  244 
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FIGURE  4 

VOCODER  ANALYSIS  OF 
BOW  ASPECT  SUBMARINE  DATA  (U) 

TW  64  W 625  Hi  5000  Hi 

SEQUENCE  14  TAPE  244 
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FIGURE  6 

VOCODER  ANALYSIS  OF 
BOW  ASPECT  SUBMARINE  DATA  (U) 


TW  = 64  W = 625  Hi  fQ  = 5000  Hi 


SEQUENCE  17  TAPE  244 
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IV.  miADKATUKK  FOKNKUTION 
(.1  . A.  {'hoot  or) 

(U  tWF'  Marly  In  thin  quarter  a computer  program  wnn  develop'd  to 

compute  a quadrat  ure  o ronscorre  I Ml  Ion , based  on  l In-  fan!  Fourier 
transform  (FKf^.  To  check  I ho  program,  a comparlnou  wnn  made  on 
Identical  not  s of  data  using  I hi'  FIT  version  and  a program  developed 
by  t ho  lUgnnl  Vt\yn  ton  Ulvlnlon.  Tho  results  woro  Identical,  and 
t hoy  ure  shown  in  Fig.  t'  of  quit  No.  ' on  Font  mot  Ni\\V4-('{l-F-  1 1 -V)  ^I'V 
Unfort  unat  oly . I ho  FF’f  voralon  wnn  embarrassingly  nlow.  A now 
version  Imn  .hint  boon  developed  that  l.akon  advantage  of  I ho  symmetrical 
proport  ton  of  tho  Four  I or  I -ransfonnat  lonn  of  real  numbers.  Tho  now 
vornlon  ( FViFl' ' In  now  opornt  lonnl , mul  Fig.  i nluavn  tho  comparative 
OvMnpnt  at  tonal  t Imon . Tho  FV’l'  progrtun  munt  oomputo  data  In  blookn  of 
poworn  of  two.  whloh  accounts  for  tho  stairstep  fund  ton  nhv'wn  In 
tho  flgnro.  Tho  two  onrvon  nlu'wtng  oomput  at  tonal  I I mo  Intornool 
tlmoo,  and  potential  nnorn  will  haw  to  rofor  to  Fig.  to  decide 
which  program  to  uno  for  eootnwy  of  tlmo.  Tho  FFI’  program  In 
written  In  lA'KTKAN  uni  does  float  Ing  point  arithmetic,  no  that  II 
can  nt  111  ho  improved . 
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FIGURE  7 

COMPUTATIONAL  TIMES  FOR  THREE  PROGRAMS 
COMPUTING  QUADRATURE  CROSSCORRELATIONS 


UNCLASSIFIED 


V.  KKIVKIV.  ANl'  rUULUjrVl’lON:'. 

(.)  . A.  ilhoot  *•  i* ^ 

(. l1  Thi-  Joint  paper  on  narrow  band  rlt-.nal  proornn l iif  l\\ 

Messra.  0.  I'.  Ora  re  and  , A.  Shooter  haa  .vet  to  l'<'  written. 

Mlaa  C.  ,t.  Webb  la  now  putting  the  flnlnhinf,  lonelier  on  a report 
dealing  with  t lie  praet leal  anpeetn  of  time  rerlen  anal.vnln  and  the 
Kl-T  a.l  Kor  1 1 hm  . 
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